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ABSTRACT 

A  calculation  procedure  is  presented  which  allows  the  one¬ 
dimensional  determination  of  flow  distributions  in  arbitrarily  connected 
(branching)  flow  passages  having  multiple  inlets  and  exits.  The  procedure 
uses  an  adaptation  of  the  finite  element  technique,  iteratively  coupled 
with  an  accurate  one-dimensional  flow  solver.  The  procedure  eliminates 
the  usual  restrictions  inherent  with  finite  element  flow  calculations. 
Unlike  existing  one-dimensional  methods,  which  require  simplifications 
to  the  flow  equations  (uncoupling  the  momentum  and  energy  equations), 
to  allow  for  arbitrary  branching  and  multiple  inlets  and  exits,  the  only 
limitation  of  the  described  methodology  is  that,  at  present,  it  can  only 
accommodate  non-rotating  configurations  (no  pumping  effects).  The 
calculation  procedure  is  robust,  and  will  always  converge  for  physically 
possible  flow.  The  procedure  is  described,  and  its  use  is  illustrated  by  an 
example. 

INTRODUCTION 

The  design  of  cooling  passages  in  turbomachinery  components 
(vanes  and  blades)  continues  to  be  a  challenging  undertaking,  since  the 
number  of  possible  cooling  configurations  is  virtually  infinite.  While  a 
number  of  capable  flow/heat-transfer  solvers  are  available  [e.g.,  ref.  1 
to  3],  they  are,  for  the  most  part,  multi-dimensional  final  design  tools  that 
require  substantial,  time  consuming,  detailed  geometric  inputs,  as  well 
as  significant  run  times.  A  need  exists  for  a  simple,  one-dimensional 
“screening”  code  that  can  be  used  to  narrow  down  the  number  of  possible 
configurations  in  a  timely  and  cost  effective  manner. 

Reference  4  describes  a  one-dimensional  computer  code  (called 
CPF,  Coolant  Passage  Flow),  which  was  developed  because  existing, 
company  pro-prietary,  one-dimensional  codes  (that  can  accommodate 
arbitrary  internal  flow  branching  and  multiple  inlets  and  exits),  make 
simplifications  to  the  governing  equations  (uncoupling  the  momentum 
and  energy  equations),  leading  to  loss  of  calculation  accuracy.  CPF  makes 
no  simplifying  assumptions  to  the  governing  one-dimensional  flow 
equations.  The  increased  calculation  accuracy,  however,  comes  at  a 


price — CPF  can  only  analyze  a  coolant  flow  path  geometry  that  consists 
of  a  single  flow  passage  with  one  inlet  and  one  exit. 

This  paper  describes  a  methodology  that  couples  CPF  with  a  finite 
element  technique  flow  solver  and  extends  CPF’s  applicability  to 
arbitrarily  connected  flow  passages  having  multiple  inlets  and  exits.  The 
main  attributes  of  CPF  are  retained,  except  for  the  ability  to  account 
for  pumping  (rotational  effects).  The  described  methodology  is,  for  now, 
applicable  only  to  cooling  geometries  in  non-rotating  components.  The 
described  methodology  is  not  limited  to  CPF,  but  can  be  implemented 
with  any  one-dimensional  flow/heat- transfer  code. 

CPF  (COOLANT  PASSAGE  FLOW) 

The  computer  code  CPF  is  described  in  ref.  4,  which  is  available  in 
the  open  literature.  CPF  was  developed  specifically  for  radial 
turbomachinery,  but  can  be  used  to  analyze  any  coolant  flow  path 
geometry  that  consists  of  a  single  flow  passage  with  one  inlet  and  one  exit. 
The  flow  path  is  defined  by  nodes  (using  Cartesian  coordinates  at  the 
passage  centroids)  and  by  intervals,  also  referred  to  as  elements.  Figures 
1  and  2  [from  ref.  4]  illustrate  the  nodal/element  representation  of  a 
typical  coolant  flow  passage  for  a  hypothetical,  cooled,  radial  turbine. 
Note  that  CPF  can  also  accommodate  tip-cap  impingement  cooling,  and 
that  a  flow  bypass  can  be  specified,  in  which  coolant  flow  is  taken  off 
at  one  point  in  the  flow  channel  and  “reintroduced”  at  a  point  further 
downstream  in  the  channel  (as  shown  in  fig.  1  and  2).  These  additional 
features,  however,  are  not  used  in  the  described,  coupled  calculation 
procedure. 

CPF  predicts  the  coolant  mass  flow  rate,  velocity,  pressure,  and 
temperature,  as  well  as  the  coolant  heat  transfer  coefficient  distribution 
along  the  passage  flow  path  from  inlet  to  exit  (for  specified  inlet  and 
exit  pressures  and  for  specified  or  empirical  friction  and  heat  transfer 
correlations).  The  solution  is  limited  to  subsonic  flow.  From  ref.  4,  the 
coupled,  one-dimensional  momentum  and  energy  equations  (without  the 
pumping  terms  and  the  mass  “deletion/reintroduction”  terms)  express 
the  rates  of  change  of  pressure  and  temperature  with  path  length  (x)  along 
an  element  as  follows: 
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Figure  1 . — Conceptual  cooled  radial  turbine  configuration 
{from  ref.  4), 
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Figure  2. — Nodal  representation  of  figure  1  upper  flow  path 
(from  ref.  4), 
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where 


A  cross  sectional  area 

Cp  specific  heat  at  constant  pressure 

Dh  hydraulic  diameter 

/  Fanning  friction  factor 

gc  gravitational  constant 

hc  heat  transfer  coefficient 

p  static  pressure 

P  perimeter 

R  universal  gas  constant 

T  temperature 

Tw  wall  temperature 

Taw  adiabatic  wall  temperature 
A'  path  length 

w  mass  flow  rate 

CPF  integrates  these  coupled  equations  (1  and  2)  along  the  flow  path. 
The  above  equations  account  for  friction,  area  change,  and  heat  transfer. 
Friction  factors  and  heat  transfer  coefficients  for  each  element  are 
obtained  from  empirically  derived  correlations. 

THE  FINITE  ELEMENT  METHOD 

The  finite  element  method  for  fluid  flow  calculations  is  described  in 
ref.  5.  Flow  passages  are  divided  into  elements  that  are  connected  at 
junction  points  called  nodes.  A  simple  stick  diagram  flow  network 
(for  this  paper’s  example  problem)  is  shown  in  fig.  3,  overlaid  on  a 
z-r  Cartesian  grid.  The  numbers  in  square  boxes  designate  the  elements, 
while  those  within  parentheses  designate  the  node  points  (shown  as  solid 
circles).  Note  that,  for  this  example,  the  number  of  nodes  (eight)  is  equal 
to  the  number  of  elements,  which  is  usually  not  the  case.  Also  note  that 
any  number  of  elements  can  be  connected  to  a  single  node.  Flow 
directions  are  assumed  (arrows  next  to  element  numbers)  in  each  element 
based  upon  known  pressures  at  the  entrance  and  exit  nodes.  In  fig.  3, 
nodes  1 ,  2,  and  5,  are  entrance  nodes,  and  nodes  7  and  8  are  exit  nodes. 
Nodes  3, 4,  and  6  are  internal  nodes,  where  the  pressures  are  unknown 
and  must  be  predicted.  The  diagram  in  fig.  3  is  an  arbitrary  connection  of 
flow  passages  and  does  not  represent  any  real  or  practical  network. 
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Figure  3. — Example  problem  stick  diagram  flow  network. 

In  the  finite  element  flow  calculation  method,  the  volume  flow 
rate  (Q)  for  each  element  is  expressed  as  the  product  of  a  “stiffness 
coefficient”,  or  “flow  function”,  (K)  times  the  change  in  pressure 
across  the  element. 

Q  =  -K-Ap  (3) 

where 

K  -  (jc»D4)/(£*h*L)  (4) 

and 

D  element  equivalent  pipe  diameter 

k  a  constant  (dependent  on  the  physical  units  being  used) 

L  element  length 

A p  pressure  drop  across  the  element 

p  fluid  dynamic  viscosity 

The  above  expression  for  Q  forces  serious  limitations  on  the  types 
of  allowable  flow  (discussed  below),  but  results  in  a  fixed  value  of  K 
for  each  element.  The  individual  FC s  are  incorporated  into  an  overall, 
global  “stiffness  matrix”,  which  is  solved  (by  matrix  techniques)  in  a 
single  step  to  give  volume  flow  rates  for  all  elements,  as  well  as  the 
previously  unknown  pressures  at  the  internal  node  points.  If  the  assumed 
flow  directions  are  correct,  the  calculated  volume  flows  are  positive.  If 
the  flow  direction  in  an  element  is  opposite  to  the  initial  guess,  the 
calculated  volume  flow  is  a  negative  number. 

At  first  glance,  the  finite  element  method  would  appear  to  be  an 
unlikely  candidate  for  flow  and  heat  transfer  calculations,  since,  in  its 
most  fundamental  form,  it  suffers  from  the  following,  serious  limitations: 
(1)  the  fluid  must  be  incompressible;  (2)  the  fluid  flow  must  be  laminar; 

(3)  each  passage  element  must  be  of  constant  cross  sectional  area; 

(4)  no  allowance  is  made  for  temperature  variations  (heat  transfer); 
and  (5)  no  differentiation  is  made  between  total  and  static  pressures. 
However,  as  will  be  shown,  all  these  limitations  are  bypassed  by  incor¬ 
porating  flow  information  generated  by  CPF.  Only  the  finite  element 
solution  method  is  utilized  in  the  described,  coupled,  iterative  solution 
procedure. 


COUPLING  OF  CPF  AND  THE  FINITE  ELEMENT  METHOD 

To  couple  CPF  and  finite  element  flow  calculations,  it  is  assumed 
that  the  “real  world”  mass  flows  through  each  element  of  a  multiple  flow 
path  can  also  be  expressed  (as  in  the  finite  element  method)  by  the  product 
of  a  “flow  function”  and  the  pressure  drop  across  the  element.  The  “real 
world”  flow  functions  must  now,  however,  encompass  all  possible  types 
of  flow  (incompressible  or  compressible,  laminar,  or  turbulent),  as  well 
as  allow  for  varying  flow  area,  friction,  and  heat  transfer  within  the 
elements.  The  values  of  the  individual  flow  functions  are,  of  course, 
unknown,  and  are  obtained  (using  CPF)  in  the  following,  iterative 
manner,  as  diagrammed  in  fig.  4. 

From  the  flow  passage  diagram  (e.g.,  fig.  3)  and  the  known  entrance 
and  exit  pressures,  flow  directions  are  assumed  (arrows  next  to  the 
element  numbers  in  fig.  3).  Next,  an  initial  guess  is  made  for  the  unknown 
node  pressures  at  the  internal  node  points,  consistent  with  the  assumed 
flow  directions  and  the  specified  entrance  and  exit  pressures.  Then,  for 
each  element,  CPF  is  used  to  calculate  the  flow  rate  and  the  “flow 
function”  (K)  by  dividing  the  calculated  flow  rate  by  the  initially  assumed 
pressure  drop.  Finally,  the  calculated  flow  functions  are  input  into  the 
finite  element  flow  solver,  which  predicts  updated  pressures  for  all  the 
internal  nodes. 

These  updated  internal  node  pressures  are  then  again  used  by  CPF 
to  calculate  further  updated  flow  rates  and  flow  functions  for  use  by  the 
finite  element  method  to  obtain  the  next  set  of  updated  internal  node 


Figure  4— Iterative  flow  calculation  diagram. 
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pressures.  The  procedure  is  repeated  until  convergence  is  achieved 
within  a  specified  tolerance  on  successively  calculated  pressures  at  all 
internal  nodes.  At  that  point  the  sum  of  all  inlet  flows  will  equal  the  sum 
of  all  exit  flows. 

As  stated  previously,  in  CPF,  empirically  derived  correlations  are 
specified  for  friction  factor  and  heat  transfer  coefficients.  These  correla¬ 
tions  can  be  specified  in  a  number  of  ways,  including  from  built-in 
correlations  or  from  curves  input  in  tabular  form.  This  feature,  (specifying 
friction  factor  and  heat  transfer  coefficients  in  a  multitude  of  ways)  is 
not  yet  available  in  the  prescribed,  coupled  solution  procedure.  Friction 
factor  and  heat  transfer  coefficients  must,  for  now,  be  specified  directly 
by  the  user  for  each  interval. 

OTHER  COUPLING  CONSIDERATIONS 

The  preceding,  simple  description  of  the  coupled  solution  proce¬ 
dure  does  not  address  the  myriad  of  details  that  are  necessary  to  assure 
a  correct  and  meaningful  flow  solution.  Only  a  few  of  the  other  consid¬ 
erations  can  be  addressed  here  due  to  paper  length  limitations.  Some  of 
the  most  important  are  as  follows: 

Pressure,  As  stated  previously,  in  the  finite  element  flow  calculation 
method,  no  differentiation  is  made  between  total  and  static  pressures.  The 
procedure,  however,  requires  decreasing  pressure  along  all  element 
flow  paths.  Therefore,  all  pressure  information  passed  between  CPF  and 
the  finite  element  flow  solver  must  be  based  on  total  conditions.  While 
static  pressure  varies  significantly  with  flow  area,  total  pressure  remains 
constant,  unless  varied  by  friction  and/or  heat  transfer.  For  the  types  of 
fluid  flow  problems  encountered  by  turbomachinery  designers,  friction 
will  always  play  a  significant  role  in  reducing  total  pressure  along  the 
flow  path. 

In  CPF,  total  pressure  is  specified  at  the  flow  path  inlet,  and  the 
corresponding  inlet  static  conditions  are  determined  from  the  calculated 
inlet  velocity.  Accordingly,  for  the  described  coupled  solution  method, 
the  known  supply  pressures  at  inlet  nodes  (entrances  to  elements),  as  well 
as  the  initial  guesses  for  the  unknown  pressures  at  the  internal  nodes  are 
specified  as  total  pressures.  The  known  exit  pressures  are  specified  as 
static  pressures,  but  only  the  corresponding,  calculated  total  pressures  at 
the  exit  nodes  are  passed  on  to  the  finite  element  flow  solver.  Also,  for 
each  element,  an  inlet  total  pressure  loss  coefficient  is  specified. 

The  finite  element  method  requirement  for  decreasing  pressure  in 
each  element  is  the  reason  why  the  described,  coupled  solution  procedure 
does  not  work  in  a  rotating  environment.  Pumping  (caused  by  rotation, 
coupled  with  increasing  path  radius)  increases  total  pressure  along  the 


passage,  and  can  easily  overcome  the  total  pressure  drop  due  to  friction. 
Efforts  to  overcome  this  limitation  have  so  far  involved  modifying  the 
stiffness  eoefficient(s)  in  the  appropriate  element(s).  While,  at  times, 
this  has  allowed  some  additional  iteration  steps  to  be  completed,  no 
overall  converged  solution  procedure  has  yet  been  achieved  in  the 
presence  of  pumping. 

Temperature.  In  the  described  coupled  solution  method,  total  tem¬ 
peratures  are  also  specified  at  all  inlets  and  as  Initial  guesses  at  the  internal 
nodes.  If  two  or  more  elements  feed  flow  into  an  internal  node,  a  “mass 
averaged”  node  total  temperature  is  calculated.  This  average  total 
temperature  is  then  used  as  the  inlet  total  temperature  for  those  elements 
flowing  out  of  the  internal  node. 

Row  Direction ,  As  stated  previously,  if,  in  the  finite  element  method, 
the  initial  guess  of  flow  direction  in  any  element  is  wrong,  the  predicted 
volume  flow  rate  will  be  a  negative  number  (although  still  correct), 
indicating  that  the  flow  is  opposite  to  the  initially  assumed  direction. 
This  is  not  the  case  with  the  coupled  solution  method.  If  an  initially 
assumed  flow  direction  is  wrong,  CPF  will  fail  to  converge  for  that 
particular  element,  since,  without  pumping,  the  flow  is  only  allowed  to 
travel  in  the  direction  of  decreasing  total  pressure.  If  this  happens, 
appropriate  error  messages  identify  the  element,  and  the  user  must  change 
the  input  to  incorporate  the  opposite  element  flow  direction. 

SAMPLE  CALCULATION 

A  sample  calculation  was  performed,  using  English  units,  for  a  * 
flow  network  based  on  the  stick  diagram  of  fig.  3.  For  simplicity,  the 
example  was  run  with  fixed  physical  gas  properties  (gamma  = 

1 .4,  Cp  =  0,24  BTU/lbm//?,  viscosity  =1.21  X 10"5  Ibm/ft/s,  gas  constant 
=  53.35  ft*lbf/lhm//?).  Zero  heat  transfer  and  zero  inlet  losses  were 
assumed,  and  the  Fanning  friction  factors  were  0.20  for  all  elements. 
Nodes  1,  2,  and  5  are  inlets,  nodes  7  and  8  are  outlets,  and  nodes  3, 4, 
and  6  are  internal  nodes  w?here  the  unknown  pressures  have  to  be 
determined.  Elements  1  and  3  have  constant  cross  sectional  area,  while 
all  other  elements  have  varying  cross  sectional  areas.  Inlet  total  pressures 
at  nodes  1,  2,  and  5  were  specified  as  50.0  psia  each,  while  inlet  total 
temperatures  at  those  nodes  were  specified  as  500.0,  400,0,  and 
300,0  °F,  respectively.  Exit  static  pressures  at  nodes  7  and  8  were 
specified  as  20.0  and  15,0  psia,  respectively.  The  initial  guesses  for 
total  pressures  at  nodes  3,  4,  and  6  were  47.0,  46,0,  and  45,0  psia, 
respectively.  Table  1  shows  the  pertinent  input  information  for  each 
element  in  tabular  form. 


TABLE  L— EXAMFLE  PROBLEM  ELEMENT  TABULATION  OF  PRESCRIBED  PARAMETERS 


Element, 

GO 

Inlet 

node, 

(n) 

Outlet 

node, 

(n) 

Inlet 

Z, 

in. 

Outlet 

2, 

in. 

Inlet 

r, 

in. 

Outlet 

r. 

In. 

Inlet 

area, 

in.2 

Outlet 

area, 

.  't 

in." 

Inlet 

total 

pressure, 

psia 

Outlet 

static 

pressure, 

psia 

Inlet 

total 

temper¬ 

ature, 

°F 

1 

1 

3 

0.0 

0,0 

-25.0 

0,0 

2.0 

2.0 

50,0 

500,0 

2 

2 

3 

-40.0 

0.0 

0,0 

0.0 

2.0 

3.0 

50.0 

400.0 

3 

3 

6 

0.0 

50.0 

0.0 

0,0 

3.0 

3,0 

4 

3 

4 

0.0 

30.0 

0,0 

50.0 

2.0 

1.0 

5 

5 

4 

30.0 

30.0 

70,0 

50.0 

2.0 

1,0 

50.0 

300.0 

6 

4 

6 

30.0 

50.0 

50.0 

0,0 

2.0 

1.0 

7 

6 

7 

50.0 

90.0 

0.0 

0.0 

2.0 

3.0 

70  O 

8 

6 

8 

50.0 

50.0 

0.0 

-50.0 

2.0 

1.0 

iv.U 

15.0 

— 
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TABLE  2.— CALCULATED  ELEMENT  PARAMETERS 
_  (Prescribed  parameters  in  parentheses) _ 


Element, 

0 

Flow 

rate, 

lbm/s 

Inlet 

total 

pressure, 

psia 

Outlet 

total 

pressure, 

psia 

Inlet 

static 

pressure, 

psia 

Outlet 

static 

pressure, 

psia 

Inlet 

total 

temper¬ 

ature, 

°F 

Outlet 

total 

temper¬ 

ature, 

°F 

Inlet 

static 

temper¬ 

ature, 

°F 

Outlet 

static 

temper¬ 

ature, 

°F 

Inlet, 

M 

Outlet, 

M 

1 

0.2349 

(50.0) 

47.14 

49.78 

46.91 

(500.0) 

500.0 

498.8 

498.6 

0.0795 

0.0844 

2 

0.2503 

(50.0) 

47.14 

49.78 

47.04 

(400.0) 

400.0 

398.9 

399.5 

0.0802 

0.0566 

3 

0.4771 

47.144 

37.56 

46.74 

37.05 

448.4 

448.4 

446.2 

444.8 

0.1115 

0.1405 

4 

0.0068 

47.144 

47.13 

47.14 

47.13 

448.4 

448.4 

448.4 

448.4 

0.0024 

0.0048 

5 

0.1884 

(50.0) 

47.13 

49.89 

46.65 

(300.0) 

300.0 

299.5 

297.8 

0.0566 

0.1210 

6 

0.1946 

47.130 

37.56 

47.00 

36.91 

305.2 

305.2 

304.6 

301.4 

0.0623 

0.1583 

7 

0.4567 

37.560 

20.84 

36.53 

(20.0) 

406.9 

406.9 

400.1 

396.8 

0.1995 

0.2425 

8 

0.2163 

37.560 

17.26 

37.34 

(15.0) 

406.9 

405.5 

405.4 

371.4 

0.0927 

0.4526 

Figure  5 —Calculated  flow  rates  (Ibm/s). 


RESULTS 

The  solution  procedure  converged  in  10  overall  iterations,  requiring 
less  than  10  seconds  of  computer  time,  using  a  550  MHz  personal 
computer.  The  calculated  flow  rates  are  shown  in  fig.  5,  while  table  2 
tabulates  further  calculated  flow  information  for  each  element.  From 
fig.  5  it  can  be  seen  that  the  sums  of  the  calculated  flow  rates  going 
into  and  out  of  the  internal  nodes  (3,  4,  and  6),  as  well  as  the  overall 
inlet  and  exit  flows  (0.6736  and  0.6730  lbm/s,  respectively),  agree  within 
the  3rd  significant  figure.  This  is  reassuring,  since  the  convergence 
criterion  is  not  based  on  flow  rates,  but  on  the  relative  differences  between 
successively  predicted  pressures  at  the  interior  node  points.  The  con¬ 
vergence  criteria  is  that,  for  each  internal  node,  ABS(PTOLD-PTNEW)/ 
PTOLD  is  less  than  or  equal  to  5X1CT4,  where  PTOLD  and  PTNEW 
are  the  previous  and  latest  calculated  node  total  pressure,  respectively. 
The  calculation  procedure  also  converged  for  other  initial  pressure 
guesses  at  the  internal  node  points  (as  long  as  they  did  not  violate  the 
assumed  flow  direction  criteria),  showing  that  the  overall  iteration 
procedure  is  robust. 


CONCLUSIONS 

1 .  An  iterative  calculation  procedure  has  been  developed  that  couples 
the  finite  element  flow  solution  method  with  an  accurate  one¬ 
dimensional  flow  and  heat  transfer  solver  called  CPF. 

2.  The  described  procedure  is  shown  to  be  capable  of  predicting  the 
one-dimensional  fluid  flow  distribution  in  arbitrarily  connected 
passages  having  multiple  inlets  and  exits  by  simultaneously  solving 
the  coupled  one-dimensional  momentum  and  energy  equations.  The 
procedure  differs  from  other,  known  methods,  which  require  simpli¬ 
fications  (uncoupling)  of  the  governing  equations. 

3.  The  calculation  procedure  can  be  adapted  to  any  one-dimensional 
flow  solver,  and  can  accommodate  any  fluid  flow  that  is  consistent 
with  the  one-dimensional  flow  code  (e.g.  compressible,  incompress¬ 
ible,  laminar  or  turbulent  flow  with  varying  flow  area,  friction,  and 
heat  transfer). 

4.  The  solution  procedure  is  robust,  and  has  always  converged  rapidly 
for  physically  possible  flow. 

5.  At  present,  the  only  limitation  of  the  calculation  technique  is  the 
inability  to  accommodate  pumping  effects.  Efforts  to  overcome  this 
limitation  are  in  progress. 
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